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Abstract 

The concept of Wardrop equilibrium plays an important role in con¬ 
gested traffic problems since its introduction in the early 50’s. As shown 
in [2], when we work in two-dimensional cartesian and increasingly dense 
networks, passing to the limit by T-convergence, we obtain continuous 
minimization problems posed on measures on curves. Here we study the 
case of general networks in which become very dense. We use the no¬ 
tion of generalized curves and extend the results of the cartesian model. 

Keywords: traffic congestion, Wardrop equilibrium, F-convergence, gener¬ 
alized curves. 


1 Introduction 

Modeling congested traffic is a field of research that has developed especially 
since the early 50’s and the introduction of Wardrop equilibrium [T^]. Its pop¬ 
ularity is due to many applications to road traffic and more recently to commu¬ 
nication networks. In our finite networks model, we represent the congestion 
effects by the fact that the traveling time of each arc is a nondecreasing function 
of the flow on this arc. The concept of Wardrop equilibrium simply says that 
all used roads between two given points have the same cost and this cost is min¬ 
imal. So we assume a rational behavior by users. A Wardrop equilibrium is a 
flow configuration that satisfies mass conservation conditions and positivity con¬ 
straints. A few years after Wardrop defined his equilibrium notion, Beckmann, 
McGuire and Winsten [3] observed that Wardrop equilibrium can be formulated 
in terms of solutions of a convex optimization problem. However this variational 
characterization uses the whole path flow configuration. It becomes very costly 
when working in increasingly dense networks. We may often prefer to study the 
dual problem which is less untractable. But finding an optimal solution remains 
a hard problem because of the presence of a nonsmooth and nonlocal term. As 
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we study a sequence of discrete networks increasingly dense, it is natural to ask 
under what conditions we can pass to a continuous limit which would simplify 
the problem. 

The purpose of this paper is to rigorously justify passing to the limit thanks 
to the theory of T-convergence and then to find a continuous analogue of Wardrop 
equilibrium. We will strongly rely on two articles [S] and [5] . The first establishes 
some hrst results on traffic congestion. The second studies the case of a two- 
dimensional cartesian grid with small arc length e. Here we will consider general 
networks in with small arc length of order e. It is a substantial improve¬ 
ment. We will show the T-convergence of the functionals in the dual problem as 
£ goes to 0. We will obtain an optimization problem over metrics variables. The 
proof of the T-convergence is constructed in the same manner as in But two 
major difficulties here appear. Indeed in the case of the grid in there are 
only four possible constant directions (that are ((1, 0), (0,1), (—1,0), (0, —I))) 
so that for all speed z G R^, there exists an unique decomposition of z in the 
family of these directions, with positive coefficients. In the general case, direc¬ 
tions are not necessarily constant and we have no uniqueness of the decompo¬ 
sition. To understand how to overcome these obstacles, we can first look at the 
case of regular hexagonal networks. There are six constant directions (that are 
exp(i{'KjQ -I- kir/S)), fc = 0,..., 5) but we lose uniqueness. Then, we can study 
the case of a two-dimensional general network in which directions can vary and 
arcs lengths are not constant. The generalization from R^ to R'^ (where d is 
any integer > 2) is simpler. Of course, it is necessary to make some structural 
assumptions on the networks to have the T-convergence. These hypotheses are 
satisfied for instance in the cases of the isotropic model in [8] and the cartesian 
one in [2]. 

The limit problem (in the general case) is the dual of a continuous prob¬ 
lem posed on a set of probability measures over generalized curves of the form 
(cr, p) where cr is a path and p is a positive decomposition of a in the family 
of the directions. This takes the anisotropy of the network into account. We 
will then remark that we can define a continuous Wardrop equilibrium through 
the optimality conditions for the continuous model. To establish that the limit 
problem, has solutions we work with a relaxation of this problem through the 
Young’s measures and we extend results in [8]. Indeed we cannot directly ap¬ 
ply Prokhorov’s theorem to the set of generalized curves since weak-L^ is not 
contained in a Polish space. First we are interested in the short-term problem, 
that is, we have a transport plan that gives the amount of mass sent from each 
source to each destination. We may then generalize these results to the long¬ 
term variant in which only the marginals (that are the distributions of supply 
and demand) are known. This case is interesting since as developed in [TJ [HI [3] , 
it amounts to solve a degenerate elliptic PDF. But it will not be developed here. 

The plan of the paper is as follows: Section 2 is devoted to a preliminary 
description of the discrete model with its notations, definition of Wardrop equi¬ 
librium and its variational characterization. In Section 3, we explain the as¬ 
sumptions made and we identify the limit functional. We then state the T- 
convergence result. The proof is given in Section 4. Then, in Section 5, we 


2 


formulate the optimality conditions for the limit problem that lead to a contin¬ 
uous Wardrop equilibrium. Finally, in Section 6, we adapt the previous results 
to the long-term problem. 


2 The discrete model 

2.1 Notations and definition of Wardrop equilibria 

Let d G N, d > 2 and fl a bounded domain of with a smooth boundary 
and £ > 0. We consider a sequence of discrete networks ilg = whose 

characteristic length is e, where is the set of nodes in fig and the (finite) 
set of pairs {x, e) with x G iV® and e G such that the segment [x, x + 
e] is included in fl. We will simply identify arcs to pairs (x, e). We impose 
\E‘^\ = max{|e|, there exists x such that (x, e) G E’^} = e. We may assume 
that two arcs can not cross. The orientation is important since the arcs {x, e) 
and (x+e, —e) really represent two distinct arcs. Now let us give some definitions 
and notations. 

Traveling times and congestion: We denote the mass commuting on arc 
{x,e) by rrf(x,e) and the traveling time of arc {x,e) by E{x,e). We represent 
congestion by the following relation between traveling time and mass for every 
arc {x, e): 

f{x,e)=g^{x,e,rrf{x,e)) ( 2 . 1 ) 

where for every e, is a given positive function that depends on the arc itself 
but also on the mass m®(a:, e) that commutes on the arc {x, e) in a nondecreasing 
way: this is congestion. We will denote the set of all arc-masses rrf{x, e) by m^. 
Orientation of networks here is essential: considering two neighboring nodes x 
and x' with {x,x' — x) and {x',x — x') G E^, the time to go from x to x' only 
depends on the mass m^{x, x' — x) that uses the arc {x, x' — x) whereas the time 
to go from x' to x only depends on the mass rrE{x', x — x'). 

Transport plan: A transport plan is a given function 7 ^ : x 1 —>■ 1R+. 

That is a collection of nonnegative masses where for each pair (x, y) G x 
(viewed as a source/destination pair), 7 ^( 0 ;,?/) is the mass that has to be sent 
from the source x to the target y. 

Paths: A path is a hnite collection of successive nodes. We therefore rep¬ 
resent a path a by writing a = (a;o, ■ • ■, with a{k) = Xk € W and 

((T(fc), a{k + 1) — cr{k)) G A® for fc = 0,..., L{a) — 1. (t( 0) is the origin of a and 
a{L{a)) is the terminal point of a. The length of cr is 


L(a)-l 

\xk+i-xk\. 

k=0 

We say that (x, e) C ct if there exists A: G {1,..., L{a) — 1} such that a{k) = x 
and e = a{k + 1) — cr(fc). Since the time to travel on each arc is positive, we can 
impose a has no loop. We will denote the (finite) set of loop-free paths by C®, 
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that may be partitioned as 


C' = u 

{x,y)GN^ xN^ 

where C| y is the set of loop-free paths starting from the origin x and stopping 
at the terminal point y. The mass commuting on the path a € will be 
denoted The collection of all path-masses w‘^(a) will be denoted w'^. 

Given arc-masses m^, the travel time of a path a G is given by: 

'rm=(o')= X! 9Hx,e,m^x,e)). 

(x,e)C(7' 


Equilibria: In short, in this model, the data are the masses 'Y‘^(x,y) and 
the congestion functions 5 ®. The unknowns are the arc-masses m^(x,e) and 
path-masses w^(a). We wish to define some equilibrium requirements on these 
unknowns. First, they should be nonnegative. Moreover, we have the following 
conditions that relate arc-masses, path-masses and the data 7 ^ : 

j’^(x,y) := w®(cr), V(a;,y) G X IV® (2.2) 


and 

m^(x,e)= ^ w^(a),V(x,e) G (2.3) 

crGC®;(ai,e)C»7 

Both express mass conservation. We finally require that only the shortest paths 
(taking into account the congestion created by arc and path-masses) should 
actually be used. This is the concept of Wardrop equilibrium that is defined 
precisely as follows: 

Definition 2.1. A Wardrop equilibrium is a configuration of nonnegative arc- 
masses : {x,e) —>■ {m^{x,e)) and of nonnegative path-masses w® : a —>■ 
w^{a), satisfying the mass conservation conditions (|2.2I) and (ED and such 
that for every {x, y) G iV® x iV® and every a G y, if w^{a) > 0 then 

< T^s(cr'),Vcr' G C^^y. 


2.2 Variational characterizations of equilibria 

Soon after the work of Wardrop, Beckmann, McGuire and Winsten [3] dis¬ 
covered that Wardrop equilibria can be obtained as minimizers of a convex 
optimization problem: 

Theorem 2.1. A flow configuration (w®,m®) is a Wardrop equilibrium if and 
only if it minimizes 

_ pm 

G‘^{x,e,m'^{x,e)) where G^{x,e,m) := / g‘^{x,e,a)da (2.4) 
subject to nonnegativity constraints and the mass conservation conditions (|2.2I) - 

ED- 
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Since the functions g'^ are nondecreasing with respect to the last variable, the 
problem (12.41) is convex so we can easily obtain existence results and numerical 
schemes. Unfortunately, this problem becomes quickly costly whenever the 
network is very dense, since it requires to enumerate all paths flows w^{a). 
For this reason, we can not use this characterization to study realistic congested 
networks. An alternative consists in working with the dual formulation which 
is 

^ -f^x,y)T^e{x,y), (2.5) 

+ {x,e)GE^ {x,y)GN^^ 

where should be understood as = {t^{x,e))(^x,e)^E^^ H^{x,e,.) := 

(G"^(a;, e, .))* is the Legendre transform of G®(x, e,.) that is 

H‘^{x, e, t) := sup{mt — G®(x, e, m)}, Vt G M+ 

m>0 

and T^e is the minimal length functional: 

Tte (x, y) = ^mn ^ (x, e). 

{x,e)(Zo' 

The complexity of (12.51) seems better since we only have = 0{e~‘^) nodes 
variables. However an important disadvantage appears in the dual formulation. 
The term T^e is nonsmooth, nonlocal and we might have difficulties to optimize 
that. Nevertheless we will see that we may pass to a continuous limit which will 
simplify the structure since we can then use the Hamilton-Jacobi theory. 


3 The F-convergence result 

3.1 Assumptions 

We obviously have to make some structural assumptions on the e-dependence of 
the networks and the data to be able to pass to a continuous limit in the Wardrop 
equilibrium problem. To understand all these assumptions, we will illustrate 
with some examples. Here, we will consider the cases of regular decomposition 
(cartesian, triangular and hexagonal, see [Figure 1[ ) for d = 2. In these models, 
all the arcs in have the same length that is e. We will introduce some 
notations and to refer to a specific example, we will simply add the letters c (for 
the cartesian case), t (for the triangular one) and h (for the hexagonal one). 
The first assumption concerns the length of the arcs in the networks. 

Assumption 1. There exists a constant G > 0 such that for every e > 0, (x, e) G 
U®, we have Ce < |e| < e. 

More generally we denote by G a generic constant that does not depend on 
the scale parameter e. 

The following assumption is on the discrete network Hg, e > 0. Roughly 
speaking, the arcs of define a bounded polyhedron (still denoted by He by 
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Figure 1: An example of domain in 2fi-hexagonal model 


abuse of notations) that is an approximation of 17. ilg is the union of cells - 
or polytopes - (Vf). Each Vf is itself the union of subcells - or facets - 
More precisely, we have: 

Assumption 2. Up to a subsequence, fig C fig' for e > s' > 0 and 17 = 
U>o^e- e > 0, we have 

17e = (^ Vf, P is finite. 

iei^ 


There exists S' G N and Si < S such that for every e > 0 and i € P , Vf = 
Convixl^^,... where for j = 1 ,... ,Sj — 1, is a neighbor of xf and 

xlis a neighbor of x^ in fig. Its interior Vfi contains no arc G . For 
i ^ j, Vf nVf = 9 or exactly a facet (of dimension < d — 1). Let us denote X( 
the isobarycenter of all nodes xf contained in Vfi. For i G P , we have 

yfi = Ffj with If finite. 


For j G If, Ffj = Conv{Xf... ,xfj^) with these {d + 1) points that are 
affinely independent. For every k I, Ff^fiFfi is a facet of dimension < d — 1, 
containing Xf. There exists a constant C > 0 independent of e such that for 
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every i € P and j G If, the volume of Vf and Ffj satisfy 

< IF/1 < Ce'^. 

The estimate on the volumes implies another estimate: for every i G P and 
fc = 1 ,..., Si, we have 

< dist(Xf, xlk) < Ce (3.1) 

(the constant C is not necessarily the same). This hypothesis will in particular 
allow us to make possible a discretization of Morrey’s theorem and then to prove 
ILemma 4.21 It is a non trivial extension of what is happening in dimension 2. 
[Figure 2] shows an illustration of [Assumption 2 jin the cartesian case. 



Figure 2: An illustration of [Assumption 2 [ in the cartesian case for d = 2. 

We must also impose some technical assumptions on and E‘^. 

Assumption 3. There exists N G N, D = {vk}k=i,...,N G and 

{ck}k=i,...,N G with a > d/p such that weakly converges in the 

sense that 

lim 1 a;, 7 ^- ) = [ ip(x,v) 0(dx,dv),yip G C(n x 

E^o+, V \e\J 
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where 9 G x is of the form 

N 

6{dx,dv) = y^Ck{x)6yf^(j.)dx. 

k^l 

The are the possible directions in the continuous model. We have to 
keep in mind that for every x G the Vk{x) are not necessarily pairwise 
distinct. The requirement Vk G with a > d/p is technical and will in 

particular be useful to prove iLemma 4.^ In our examples, the sets of directions 
are constant and we have 


Dc = {uci = {1,0), Vc2 = (0,1),?^C3 = (-l,0),?;c4 = (0,-1)}, 
Dt = Dh = {vt, = e*"/® • 


The Cfc’s are the volume coefficients. In our examples, they are constant and 
do not depend on fc: 

c, = l.c, = ^a„dc, = ^. 

We notice that the c/s are different. Indeed, a square whose side length is e 
does not have the same area as a hexagon whose side length is e. The next 
assumption imposes another condition on the directions VkS. 

Assumption 4. There exists a constant C > 0 such that for every (x,z,^) G 
X X there exists Z G R.+ such that \Z\ < C and 


N 

Z ■ ^ = min{Z ■ Z = {zi,..., zn) G and ZkVk{x) = z}. (3.2) 

fe=i 


This means that the family {u^} is positively generating and that for every 
{x,z) € R'^x G R'^, a conical decomposition of z, not too large compared to z, 
is always possible in the family D{x). In the cartesian case, for all z G R^, we 
have an unique conical decomposition of z in Dc while in the other examples, we 
always have the existence but not the uniqueness. The existence of a controlled 
minimizer allows us to keep some control over ’’the” decomposition of z G R"’* 
in the family of directions. We now see a counterexample that looks like the 
cartesian case. We take N = A, d = 2, vi = (1,0),= (—1,0) and = (0,1) 
(these directions are constant). We assume that there exists Xq in such that 
V 2 {x) —>■ —V 4 = (0,1) as x —>■ xo and U2i(a::) > 0 ior x ^ xq where ^2(2;) = 
(f2i(a;),'y22(2^))- Then for every z = (zi,Z2) G R^ such that Zi > 0, we can 
write z = X 2 {x)v 2 {x) + X 4 {x)v 4 {x) for x close enough to xq with 


A2(a;) 


Zl 

V2l{x) 


and A4(a;) = zi 


V22{x) 

V2l{x) 


- Z2- 


Then A2(a::) and A4(a;) —>■ +00 as x —>■ xq. For ^ = (1,0,1,0), the value of 
always is zero but the only decomposition that solves the problem is not 
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controlled. This example shows that the existence of a controlled decomposition 
does not imply that the minimal decomposition is controlled. However this 
assumption is still natural since we want to control the flow on each arc in order 
to minimize the transport cost. 

Assumption 5. Up to a subsequence, may be partitioned as = El 
such that for every k = 1,... ,N, one has 


lim 

E^0+ 


E 

{x,e)GEl 


|e|V 



and for every (x, e) € Ef,, 


Ck{x)(p{x,Vk{x)) dx, V(^ S X ^,R). 


(3.3) 


n “ 

e 


= 0 ( 1 ). 


(3.4) 


This hypothesis is natural. Indeed, for every x € fl the condition Ck{x) > 
0 implies that for e > 0 small enough, there exists an arc {y, e) in E^ such 
that y is near x and e/|e| close to Vk(x). We will use it particularly to prove 
ILcmma 4.41 The next assumption is more technical and will specifically serve 
to apply iLemma 4.3l 

Assumption 6. For e > 0, there exists d (finite) sets of paths Cf,...,Cl 
and d linearly independent functions ei,..., H —>■ such that for every 

X G = l,...,d, ei{x) = J2k‘^k'^k{x)vk{x) where for k = 1,...,N, is 
constant and equal to 0 or 1 so that for i = 1,... ,d, we have 


[J {(a:, e) C cr} = {(a:, e) G E'^/3k G {1,..., N}, a\ = 1 and {x, e) G El}. 

(TgCf 


For every {cr,a') G Of x Of, if a ^ a' then cr fl tr' = 0. We assume 


max (dist(yQ,dn), dist{yj^(^cr)j9n)) —>■ 0 as e —>■ 0 . 


and 

\\yk+i - yk\'^~^ - \yk - yk-i\’^~^\ = 0{e^) 
for a = {yo,.. .,yL{cr)) S Of and k = 2,.. .,L{a) - 1. 

Formally speaking, it means that we partition points into a set of disjoint 
paths whose extremities tend to the boundary of H and such that we may do 
a change of variables for the derivates. For ip a regular function, then 
becomes {d^^ip ,..., de^^ip). Note that the condition on the modules is a sort of 
volume element in spherical coordinates, dV being approximately r‘^~^dr. In 
our examples, this requirement is trivial since arcs length is always equal to e in 
Hg. It will allow us to prove some statements on functions. More specifically, we 
will need to show that some functions have (Sobolev) regularity properties. For 
this purpose, it will be simpler to use these e^. In the cartesian case (d = 2), we 
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O’!,!; ■ ■ ■ j <^1,5 S C'l 

0'2,1, ■ ■ ■ , 0-2,5 £ C'l 


n 


Figure 3: An illustration of [Assumption ^in the hexagonal case for d = 2 


can simply take and Ccj = Vc^ and in the triangular one, "we can -write 

ejj = CtVti and ejj = CtVt 2 - In the hexagonal case, we can for instance define 
Cft,! = Ch{vhi + Vh^) and = Chivu^ + Vh^,)- This last example is illustrated 
by [Figure 3| (with 17 being a circle). We can notice that some arcs are in paths 
(T G Cf and cr' G C|. 

All these structural hypothesis on are satisfied in our three classical exam¬ 
ples. The cartesian one is the most obvious and the hexagonal one is a subcase 
of the triangular one. The following assumption is on the transport plan. 

Assumption 7. (e^“^7^)e>o weakly star converges to a finite nonnegative mea¬ 
sure j on X n in the sense that the family of discrete measures 
J2{x,y)&N'^xN-^ Weakly star converges to 7 ; 


lim V £2 ^j^{x,y)(p{x,y) 

e—>- 0 + 

{x,y)^N^^ 



ipd'y] 


Vyi G C(17 X 17). 


The next assumption focuses on the congestion functions g^. 
Assumption 8. g^ is of the form 

m 


5®(a;,e,TO) = \e\‘^^'^g ( x, 


Ve > 0, (x, e) G if®, m > 0 


(3.5) 


(3.6) 


where (7 : 17 x ^ x R+ 1 -^ M. is a given continuous, nonnegative function that 
is increasing in its last variable. 
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We then have 


( 6 TTl \ 

cc, —’ I 1^/2 ) ■''^here G'(a;, m) := J g{x,v,a)da 


and 


H^x,e,t) = \e\‘^H 


x - 

’ \ene\^/y 


where H{x, v, •) 


(G(x,v,-))* 


i.e. for every ^ G K+ : 


iJ 



sup 

me]R+ 


— G 



For every (x, v) G fi x H{x, v, •) is actually strictly convex since G{x, v, •) 

is G^ (thanks to [Assumption 8|). Note that in the case d = 2, we have 


g^x,e,m) 


|e|5 



for all e > 0,{x,e) G E^,m> 0. 


It is natural and means that the traveling time on an arc of length |e| is of order 
|e| and depends on the flow per unit of length i.e. m/|e|. For the general case, 
we have extended this assumption. The exponent d/2 is not very natural, it 
does not represent a physical phenomenon but it allows us to obtain the same 
relation between G and H, that means, H{x,e, •) is the Legendre transform of 
G{x,e, ). Moreover, we may approach some integrals by sums. We can also 
note that there is actually no e-dependence on the g^. 

We also add assumptions on H: 

Assumption 9. H is continuous with respect to the first two arguments and 
there exists p > d and two constants 0 < A < A such that for every (x,v,^) G 
n X X ]R_|_ one has 


x{e-i)<H{x,v,o<Ke+ 1)- ( 3 . 7 ) 


The p-growth is natural since we want to work in in the continuous limit. 
The condition p > d has a technical reason, that will allow us to use Morrey’s 
inequality. That will be crucial to pass to the limit in the nonlocal term that 
contains Tfi ■ 


3.2 The limit functional 

In view of the previous paragraph and in particular [Assumption 8[ it is natural 
to rescale the arc-times t® by defining new variables 

^^(^> (^’ ^ (^■^) 
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i.e. for every (x, e) G , 


ce, ^ S®(a^,e,TO®(x,e)) / e m‘^{x,e) 

? e) = - rUT^ - =9\x,—, 


V |er 

Then rewrite the formula (j2.5p in terms of as: 

inf , J^r) := 

eeRj 


where 


and 


(3.9) 


mi--= E |e|"i^(x,i^,e(2:,e)) (3.10) 

{x,e)^E^ \ M / 


E nun ^ |e|''/2e(^,e) . (3.11) 


, o-eCi „ 

{x,y)^N^‘^ y ’ (j2,e)Ccr 

In view of |Assumption 9[ let us denote 

LI{9) := {C G LP{n X S'^-\6»),5 > 0}. 
It is natural to introduce 


/o(e) := [ H{x, V, e(x, u))0(dx, du), G L^(d), (3.12) 

dnxS-^-i 


as the continuous limit of Iq. It is more involved to find the term that plays the 
same role as If since we must define some integrals on paths. First rearrange the 
second term (13.1111 . that is, as an integral. Let G , (x, y) G iV® x and 
a G Cf y, a = (xo,..., Xj;,(o.)). Let us extend a on [0, L(cr)] in a piecewise affine 
function by defining a{t) = a{k) + {t — k){a{k + 1) — cr(fc)) and in a piecewise 
constant function : ^^{a{t),&{t)) = ^^{a{k),a{k + 1) — cr(fc)) for t G [k,k + 1]. 
For every (x, e) G we call 'I'®(x, e) the ’’canonical” decomposition of e on 
D{x). More precisely, recalling [Assumption 5[ for (x, e) G E^, there exists 
k(x,e) G {!) ■ • ■) fV} such that (x, e) G Ef^^ and then we set 

«'®(x,e) = ( 0 ,..., |e| ,..., 0 )gR^. 

^(x,e)th coordinate 


For a £ and t G [k, k + l[, we write = 'I''^(tT(/c), CT(fc + l) —tT(/c)). 

Let us also define a function as follows: 

^^(x) = E ^-4^7-^d>®(x, e) for X G A^^. 

e/{x,e)GE^ 
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Let us extend in a piecewise constant function on the arcs of E^: let y G 
such that there exists (x, e) G with y G (x, e), then we define 


s-te) = E 


e) 


’i'^(x, e). 


(x,e)G-E^/yG(2^,e) 


This definition is consistent since the arcs that appear in the sum are in some 
different E^. By abuse of notations, we continue to write cr, and for these 
new functions. Thus we have 

X! |e|r(a;,e)= ^ lcr(k + 1) - a(k)l^^(a(k), (T(k + 1) - (T(k)) 

(x,e)C<7 k—0 

L(cy)-1 

= ^ 'I'^(a(k),cr(k + 1) - a(k)) ■ ^^(a(k)) 
k^O 

= [ ■ ^^{a{t))dt. 


We then get 

|e|f (x,e) = ^Jnf ■ ^^{ait))dt 


{x,e)(Zo 


= inf / ^^d(t),d{t)) ■ ^%d{t))dt 

'^eC^x.y Jo 

where d : [0,1] —>■ is the reparameterization of (t(t) = a{L{a)t),t G [0,1]. For 

every a; G z G let us define 

Al = Iz G IR+, Z = (zi, Zn)/ Y ^kVkix) = 2:| . 


k=l 


Then for every ^ G C{n x ^,R+), define 


c^{x,y)= inf 

O X , y 


•^0 \k=l J ) 

per, Jq J 


= jnf inf L^{(t,p), 

<7€Cx,y 


where 


= |p : t G [0, 1] -J> p{t) G IR+ /CT(t) = Y^ Vk{(j{t)) pk{t) a.e. 11 
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Li{(r,p) = Pk(t) 


dt 


and Cx,y is the set of absolutely continuous curves cr with values in Q and such 
that cr(0) = X and ct(1) = y. For every x G U and z G R'^, there is no a priori 
uniqueness of the decomposition of z in the family {vk{x)}k so that we have to 
take the infimum over all possible decompositions. The definition of Va takes 
in account this constraint. For every p and t G [0,1] the terms pk{t) are 

the weights of a{t) in the family {vk{cr{t)}. {a, p) is a sort of generalized curve. 
It will allow us to distinguish between different limiting behaviors. A simple 
example is the following approximations, e > 0 : 


ffizvw: 


It is the same curve a in both case but the p are different. Indeed, in the 
left case, we have only one direction that is (1,0), it is a line. In the right one, 
we have two directions that tend to (0,1) and (0, —1) as e tends to 0, these are 
oscillations. 

Now, our aim is to extend the definition of cj to the case where ^ is only 
L^{6). We will strongly generalize the method used in [5] to the case of gener¬ 
alized curves. First, let us notice that we may write cj in another form: 


C((x,y) = inf Lf(a) for ^ G C(fl x ^,R-i-), 


where ^ 

Zj(cr) = /" dt (3.13) 

Jo 

with for all a: G H and y G ${(a:, y) being defined as follows: 

$ 5 (x, 2 /) = inf ^yfe^(a;,Ufc(x)) : Y = {yi,--- ,2/jv) e 
U=i 

= inf y • ^(x), 

Y&Al 

where $,{x) = {^{x,vi{x)),... ,^{x,VNix))). 

The next lemma shows that $5 defines a sort of Finsler metric. It is an 
anisotropic model but <I>^ is not even and so cj is not symmetric. Moreover 
cj is not necessarily positive between two different points so that is not a 
distance. However ${(a:, •) looks like a norm that depends on the point x G ft. 
Its unit ball is a polyhedron in R^^ that changes with x. Formally speaking, the 
minimizing element Y = {yi, ..., j/jv) represents the coefficients for the ’’Finsler 
distance” c^. 

Lemma 3.1. Let ^ G C{n x S'^^^jR^). Then one has: 

1. The function is a minimum and continuous. 
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2. For all x € fl, •) is homogeneous of degree 1 and convex. 

3- If {fn} is a sequence in C{fl x and x G H such that {^n(a^)} 

converges to ^{x) then converges to ^^{x,y) for all y G R'^. 

The proof is left to the reader. 

Due to ILemma 3.11 we may reformulate (13.131) : 

^d<^it),^{t))dt = \&it)\^^(^a{t),-^^^dt. 

The following lemma gives a Holder estimate for cj and extends Proposition 3.2 

in[H]. 

Lemma 3.2. There exists a nonnegative constant C such that for every ^ G 
(^(D X and every {xi,X 2 ,yi,y 2 ) G , one has 

Icd^uVi) - Ci{x 2 ,y 2 )\ < C||^||iP(e)(|a;i - X 2 f + \yi - y 2 f), (3.14) 

where /3 = 1 — d/p. So if {fn)n G C(fd x is bounded in U'{9), then 

(c{„)n admits a subsequence that converges in C{il x r2,R+). 

For every f G L^{0), let us define 

cd^,y) = sup{c(a;,y) : c G ^(0}, V(x, y) £ D x D (3.15) 

where 

^(^) = {limc^„ in C(n x D) : (en)n £ C(D x §^-\R+)",^„ ^ ^ in LP{e)). 

n 

We will justify in the next section that for every ^ £ L?J.(6(), the continuous 
limit functional is 

J{0 ■= loiO - IiiO = [ H{x,v,f{x,v))6{dx,dv) - (_ _c^d'y. (3.16) 

The two following lemmas are generalizations of Lemmas 3.4 and 3.5 in [5]: 

Lemma 3.3. If^ G L+(0) then there exists a sequence {fn)n m C'(Dx§‘^“^,R+) 
such that converges toc^ in C{ft x fl) as n ^ oo. 

Lemma 3.4. If G C(D x S'^“^,R_|_) then cj = C{- 

3.3 The F-convergence result 

We will prove that the problem (13.161) is the continuous limit of the discrete prob¬ 
lems dSll) in the T-convergence sense. The T-convergence theory is a powerful 
tool to study the convergence of variational problems (convergence of values but 
also of minimizers) depending on a parameter. Here we want to study problems 
depending on a scale parameter (which is e), it is particularly well suited. Ref¬ 
erences for the general theory of T-theory and many applications are the books 
of Dal Maso [10] and Braides [5] . 

First let us define weak convergence of a discrete family £ R^^ . 
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Definition 3.1. For e > 0, let and ^ S L^{9), then is said to 

weakly converge to ^ in ^ if ■ 


1. There exists a constant M > 0 such that Ve > 0, one has 

liriU,p:=( E <M, 

\{x,e)eE'^ j 

2. For every ip G C{td x one has 

lim E M‘^Fix,-^]c{x,e)= <p{x,v)£,{x,v)9{dx,dv). 


{x,e)GE^ 


Definition 3.2. For e > 0, let F® : —>■ R U {+oo} and F : F+(0) —>■ 

RU {+oo}, then the family of functionals (F^)^ is said to T-converge (for the 
weak topology) to F if and only if the following two conditions are satisfied: 


1. (T-liminf inequality) \/f G L^{9),f^ such that —>■ one has 

liminfF"(e)>F(0, 

£->■ 0 + 


2. (T-limsup inequality) G Lf^{9), there exists ff G R^'® such that 
and 

limsupF''(e) <F{f), 

E— >0+ 

Now, we can state our main result, whose complete proof will be performed 
in the next section: 


Theorem 3.1. Under all previous assumptions, the family of functionals T- 

converges (for the weak topology) to the functional J defined by (13.161) . 

Classical arguments from general F-convergence theory allow us to have the 
following convergence result : 

Corollary 3.1. Under all previous assumptions, the problems (|3.9I1 for all e > 0 
and (j3.16ll admit solutions and one has: 


min 




min 


AO- 


Moreover, if for any e > Q, is the solution of the minimization problem (1331) 
then where f is the minimizer of J over L\(9). 
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Proof. First, due to (1X71) we have 

/o^(r)>A ^ |e|'^(f(x,er-l) 

{x,e)GE^ 

= Ani?,,-A^|er 

since from [Assumption 5| it follows that 

N 


E 

{x,e)GE^ 




f 6{dx,dv) ='S^ [ Ci{x)dx = C. 
JnxS''-! k=i'^^ 


max 

e/(-.e)eB' 


s,p 


To estimate the other term let us write 

d^{x,y)= nun ^ |e|^^(z,e). 

(z,e)Co- 

Vxq S N^,\/x,y G two neighboring nodes, we have 

\c^{xo,x) - P{xo,y)\ < max |e|^^(x,e) 

e/(a;,e)eE= 

<£ max £f{x,e). 

ej {x,e)^E'^ 

Then thanks to lLemma 4.21 we have for every x,y G N^, 

\d'{xo,x) - c®(a;o,y)| < C 
hence with xq = x 

c^{x,y) < C 
< 

so that recalling (13.111) . we have 

\mn\<cmkp e‘^/^-^^kx,y)<cu^\u,„ 

x,y£N^ 

because it follows from |Assumption~7| that 

f dy as £ —)> O’*". 

x,yGN^ ExQ 

In particular, we get the equi-coercivity estimate 

j-(a > - riu.p -1). 


max £^(-,e) 

e/(.,e)G£;= 
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Since J‘^ is continuous on , this proves that the infimum of over 
is attained at some and also that ||^'^||£,p is bounded, in particular, we can 
define for e > 0 the following Radon measure Mg : 


(a;,e)G-E= ^ 

Due to Holder inequality we have for every e > 0 and ip G C{^1 x 

|(M„(^)|<C||<^|U,, (3.17) 

where q = p/{p — 1) is the conjugate exponent of p and the semi-norm ||.||£,g is 
defined by: 

/ q\ 

Because of [Assumption 3| there is a nonnegative constant C such that < 

C’||‘^||oo for every p S C{Q x We deduce from (13.171) and Banach- 

Alaoglu’s theorem that there exists a subsequence (still denoted by Mg) and a 
Radon measure M over fl x with values in K to which Mg weakly star 
converges. Moreover, due to (13.171) and [Assumption 3| we have for every p G 
C{Ti X §^-i,R) 

|(M,^)| < C lim W^pWe^q = C\\p\\li{b) 

£—>•0 + 

which proves that M in fact admits an LP(0)-representative denoted by ^ G 
L\{0) (componentwise nonnegativity is stable under weak convergence) and 
^ in the sense of idefinition 3.11 We still have to prove that ^ minimizes J 
over LFj^(9). First from the F-liminf inequality we find that 

J{0 — lim inf ) = lim inf min J^(^®). 

£->0+ £->0+ JegR#®' 

Let C G L^{9), we know from the F-limsup inequality that there exists a se¬ 
quence (C^)£ G such that (j® —>■ C in the sense of idefinition 3.11 and that 

limsup J^(C) < J(C). 

£->■0+ 

Since minimizes we have that 

J(e) < liminf J'^(C^) < limsup J'^(f) < limsup J^(0 < J(C). 

e^ 0 + g_,.o+ £_>o+ 

We can then deduce that ^ minimizes J over L+(0) and we have also proved 
the existence of a minimizer to the limit problem. We also have that 

min J < liminf min < limsup min J‘^{C) < '^(C); ''^C ^ 

Li|.(e) £-> 0 + .tegK*®' £-> 0 + 5 'eR*®' 
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which provides the convergence of the values of the discrete minimization prob¬ 
lems to the value of the continuous one. Furthermore we have convergence of 
the whole family and not only of a subsequence by the uniqueness of the min- 
imizer ^ of J over L^(9) since J (and in particular Iq, due to [Assumption 8| ) is 
strictly convex. □ 


4 Proof of the theorem 


4.1 The F-liminf inequality 

For £ > 0, let G and ^ G L^{9) such that ^ (in the sense of 

Idefinition 3.11) . In this subsection, we want to prove that 

liminf J"(e) > J(0. (4.1) 

e^0+ 

We need some lemmas to establish this inequality. The first one concerns 
the terms Iq and Ip. 

Lemma 4.1. One has 

hminf/^(^'') > /o(0- 
£->0 + 

The proof is similar to that of Lemma 4.2 in [2]. Now we need the following 
discrete version of Morrey’s inequality to have information about the nonlocal 
term. 

Lemma 4.2. Let 6^ G and G such that 

\p‘^{x) — < e9^{x), for every x G and every y neighbor of x, (4.2) 


then there exists a constant C such that for every {x,y) G E‘^ x E‘^, one has 


\ip%x) - ip%y)\ < C\\ 9 ^\\e,p{\xi -yi\ + \x2- y2\f 


where /3 = 1 — d/p and 


11^11 


£,p — 



Proof. The idea is to linearly interpolate in order to have a function in 
VF^’P(n) and then to apply Morrey’s inequality. By recalling [Assumption 2[ 
let = Conv(a;f, ■■■,x‘f ) a polytope in the discrete network where for k = 
1,..., L, x% is an neighbor of x‘/_^_i in N^, the indices being taken modulo L. 
Let us denote the isobarycenter of all these nodes (it is in the interior of 
by assumption) and let us define 


= E 


L 
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Let {Fj} denote the subpolytopes given by [Assumption 2| for with Fj = 
Conv(X'^, and these {d + 1) points that are linearly independent. 

Then for every x in 1/^, there exists j such that x € F?. a; is a conical combi¬ 
nation of ... ,Xj^. There exists some unique nonnegative coefficients 

A, Ai,..., Ad > 0 such that x = XX^ + J2i=i 

d 

= A^(x-) + ^ {XI) . (4.3) 

i=l 

We still denote this interpolation by We then have G with 

||V(/5®||p < C||d'^||e_p (computing is technical and we detail it below). We 
conclude thanks to Morrey’s inequality. 

Computing Vt/?®: 

We take the notations used in the above proof but we remove the e-dependence 
for the sake of simplicity. Taking into account the construction of tp®, we will 
compute on a subpolytope F = Conv(X,Xi, ...,Xd) of V where X is the 
isobarycenter of all nodes Xk in and X, Xi ,..., Xd are linearly independent. 
We rearrange the Xk’s such that Xk is a neighbor of Xk+i- Let H be the affine 
hyperplane 


n = 



X 


Xi 

^^{Xi) 


Xd 

‘P^Xd) 


i.e. the affine subspace containing {X,(fi^{X)) and directed by the vector space 
generated by the vectors {Xi — X,ip^{Xi) — tp^{X )),..., (Xd — X,tp^{Xd) — 
(p^{X)). It is a hyperplane since the points X,Xi,...,Xd are linearly in¬ 
dependent. Then there exists some constants ao,...,ad+i such that 'H = 

{(21 ,..., Zd-i-i) G -hOd-i-i^d+i-l-ao = 0} with Od-i-i ^ 0 (otherwise 

X, Xi,..., Xd would not be independent). The normal vector to the hyperplane 
is n = (oi,..., Od+i) and 


= 


( -ai/ad+i 


V —0,d/o,d+l 


Without loss of generality we assume X = 0, Xi = (j/i, 0), X2 = (j/J, 2/2 ,0),..., Xd 
{yf,yd), with j/i,..., j/d ^ 0. Then we have for /c = 1,..., d 


( Xk-X 
{ ip%Xk) - ip^iX) 


■n = 0 = ^ -k akXk + ad+i{p'"{Xk) - </?®(X)), 
i<k 


I.e. 


Qfc - _2j_y±_ p^Xk) - p^jx) 

ad+i ^ ad+i yd yk 


(4.4) 
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We must now find an estimate on each of these terms. First for k = 1,... ,d, 
there exists ik G {1,..., L} such that Xk = Xi^ and it follows from (14.21) that 
we have 




v^ix)\ = 




1^1 

L 


< £^ 6 ^{ xi ) 


Z=1 


with the indices taken modulo L. Moreover, due to m, for k = 1,... ,d and 
i < k, we have \yf\ < |^ — < Ce. Taking into account the fact that the 

XkS are linearly independent, we get \yk\ > Ce. Then from (14.41) it follows for 
k = 1,..., d 

ak 

Od+l 




Od+l 


cY^e^xi). 


(4.5) 


We finally obtain by an induction on k that 


Qfc 

O-d+l 


L 

<cY^e%xi) 

1=1 


so that due to (14.3|) . |Vv5^(a:)| < C^2i for every x in the subpolytope 

F. We finally conclude with | Assumption~2] that 

l|v^^||p<q|riu,p, 


which completes the proof. □ 

The discretization of the Morrey inequality is crucial since we may now 
extend c® on 17 x 17. For every (x, y) G x we define 

/■i 


C{x,y)= nhn ^ |e|^''(2:, e) = inf / ^^{d{t),d{t)) ■ $,^{d{t))dt. {A.Q) 

(Z,e)ca 

By definition, if xq £ 17e and x and y neighbors in 17^, we have 
C{xq,x) < C{xQ,y) + £ max e). 

e/{y,e)eE’^ 


Since ||C®||e,p is bounded, we deduce from ILemma 4.2l that there exists a constant 
C such that for every £ > 0 we have 

\F{x,y) - c®(xo,?/o)| < C{\x-xof + \y- yof), V {x,y,xo,yo) £ 

We can then extend C to the whole 17 x 17 (we still denote by C this extension) 
by 

F{x,y):= sup {F {xo,yo) - C{\x - xo\^ + \y - yof)}, ^ ix,y) G x n. 
ixo,yo)eS^eX£le 
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By construction, still satisfy the uniform Holder estimate on the whole QxQ 
and since c® vanishes on the diagonal of H x fl, it follows from Arzela-Ascoli 
theorem that the family (c^)e is relatively compact in C{il, x H). Up to a 
subsequence, we may therefore assume that there is some c S (7(17 x U) such 
that 

c® —>■ c in (7(17 X 17) and c(x, x) = 0 for every a; £ 17. (4-7) 

From [Assumption 7| it may be concluded that 

e^-^c‘^ix,y)j‘^ix,y)-)■ [ _cd'r. 

(a:,y)eAf=xAf= dOxfi 


In consequence, with ILemma 4. 11 it remains to prove c < cj on 17 x 17. We will 
show that c is a sort of subsolution in a highly weak sense of an Hamilton-Jacobi 
equation and we will then conclude by some comparison principle. The end of 
this paragraph provides a proof of this inequality. 

Lemma 4.3. Let xq £ 17, ^ tp £ IU^’^(17) such that ipixo) = 0 

(which makes sense since p > d so that p is continuous). If for a.e. x G LI one 
has 


Vp{x) ■ u < 4)^(a;, u) 


inf 


Xkf,{x,Vk{x)) 


for all u gM.‘^ 


(4.8) 


then p < c^(xo, •) on 17. 

Remark 1. The above assumption (|4.8I) is equivalent to : 


L/p{x) ■ Vk{x) < ^k{x) = £,{x,Vk{x)), \/x a.e.yk = 1,..., W 

Proof. The result is immediate iip G (17) and ^ is continuous on 17. Indeed, in 
this case, assumption (|4.8p is true pointwise and if a; £ 17 and cr is an absolutely 
continuous curve with values in 17 connecting xq and x then by the chain rule 
we obtain 

Vp{a{t)) ■ &{t)dt < f ^eio'{t),(7{t)) dt 

Jo 

and taking the infimum in a we get p < cj(a;o, •) on 17 so that p < c^(xo, •) due 
to lLemma 3.41 For the general case, if p is only W^’P{LI) and ( only L^(6(), we 
first extend p to & function in IU^’P(K.'^) and we extend ^ outside 17 by writing 
f{x,v) = |V(^|(a;) for every x £ M.^,v £ so that if a; £ R'^\17 and u G S‘^~^ 
we have 



Vp(x) ■ u < |V(/ 3 (a:)| 

< \Vp{x)\\U\i = ${(a:,u), 

where U = (mi, ..., un) £ A“ is a minimizer of ^^(x, u) The fact U G Af. implies 
that |m| = 1 < \U\i. By homogeneity of (14.8p in u, (14.811 thus continues to hold 
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outside ri with the previous extensions. We then regularize yi and Let us take 
a mollifying sequence Pn{x) = n‘^p{nx), x G where p is a smooth nonnegative 
function supported on the unit ball and such that P = 1- Set := Pn 
and ipn '■= Pn* — {pn * Let x G Recalling [Assumption 4| (Holder 

condition on the Vks), we have 


V(p„(x) • Ufe(x) = / pn{y)y‘p{x - y) ■ Vk{x) dy 

< / Pn(y)ik(x - y) dy 

+ / Pn(2/)V(p(x - y) • (ufe(x) - Ufc(x 

jR-i 

< ^(((x)+n‘^““||p|loo / |V(p(x 

<e,”(x) + Cn'^-“-''/^||p||oo||Vp||p 

= ^k(x) + 


2/)) dy 
y)\dy 


where e„ > 0, —>■ 0 as n —>■ oo (since a > d/p). So by using the above remark 
and the previous case where (p and ^ were regular, we have < c^n_|_£^(xo,.) 
and from the convergence of to p it follows that 


p = limsuppn < limsupc^n+£^(xo, •) < c^{xo,-), 

where the last inequality is given by the definition of as a supremum (13.151) 
and the relative compactness of cj"+e„ in C{Q, x H). □ 

We want to annlv lLemma 4.31 to c(xq,-) so that we need c(xo,') G 
for every xq G H. Let (ei,...,ed) given by [Assumption 6} <p G and 

Xq G fie such that |xo — xgl < £■ Using the uniform convergence of c®(xg,-) to 
c(xo, •) and I Assumption 5} for ip G C^(H) and 2 = 1,..., d we have 


Tip := / c(xo, x)Vip(x) • ei(x) dx 

Jn 


crGCf k=0 lyk+L ykl 


,<p(yfc+i) - piyk) 
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where a = {yo ,..., yL{<j))- Then we can rearrange the sums as follows 


T^^p = 

/Her)-! 

lim^ ^ X! <f{yk) {\yk - yk-if~^c‘^{xo,yk-i) - \yk+i - ykf~^c^xo,yk)) 

crGCf \ k=l 

+ ‘P{yLia)) \yL{a) “ y C^Xq , y L{a)-l) “ ^{yo)\yi “ 2/0 I(a^O , 2/o)^ 

L(<t )-1 

X! ‘Piyk){\yk-yk-if~^c''ixQ,yk-i)-\yk+i-ykf~^c''ixQ,yk)) 

crGCf A:—1 

since for e small enough, yo and 2/L(cr) not in the support of (p thanks to 
[Assumption 6[ For a € Cf, we thus have 


{\yk - yk-if~^c^{xo,yk-i) - ly^+i - ykf~^c^{xo,yk)) 

k^l 

= X! (v^(2/fc)[l2/fc “ 2/fc-il'*~^(c''(a;o,2/fc-i) - c''(a:o,2/fc)) 

+ c^(a;o,2//c)(l2/fe - yk-if~^ - \yk+i - ykf~^)]) 

L(a)-1 

<Ce’^ \ip(yk)\- 

k^l 

Indeed in the first term, we use the fact that if x and y are neighbors in fig then 

c®(xo,a;) < c^{xQ,y) + |a: - y| max^®(y, •) 

and we obtain the upper bound on the second term due to [Assumption 6| Hence 
by using Holder and the fact that ||^®||e,p is bounded, we obtain 


/ c{xq, ■)Vlp ■ ei 




JQ. 

This proves that c(a;o,-) G By a similar argument we obtain that 

c(‘) yo) G for every yo G H. 

Lemma 4.4. Let xq G LI and c be defined by (14.7F one has 
1. For every w G the following inequality holds 


/ Vxc{xq,x) ■ w{x) dx < / ^^{x,w{x)) dx. (4-9) 

an an 

2. c < and so one has the T-liminf inequality. 
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Proof. 1.) Let a{x) = {ai{x), ...,aNix)) be a minimizing decomposition of 
w(x) i.e. for all x G 


N 


N 


XGA: 


inf 'Y' Xk^k{x) = Y] ak{x)f,k{x) 

.w(x) ^ 


k=l 


fc=l 


with of course w{x) = ^ ak{x)vk{x) and ak{x) > 0. Then we have 

/ Vxc{xo,x) ■ w{x) dx = 'Y / ak{x)Wxc{xo,x)-Vkix)dx. 

J n k—i ^ 


However the a^’s are not necessarily smooth so we must regularize the ak to 
pass to the limit. As usual, we consider a mollifying sequence {p^) (with d > 0), 
write 

N 

<^l = and 10 ^=^^ aivk 

k^l 


for fc = 1,..., A^. Hence we have 


Vxc(xo,-) 


' w = lim 
<5^0 


VxC(xo,-) • 


Let Xq G such that |xo — Xq| < e so that we have the uniform convergence 
of c®(xq,-) to c(xo,-)- By using [Assumption 5[ for every (p G we know 

that for fc = 1, N, 


Ck{x)p{x)V xc{xo, x) ■ Vk{x) dx 


= lim 

£—>•0 + 


E 

(x,e)GEl 


+ e) - P{xl,x) 


So we may write for a fixed 5 


N 


Vxc{xo,-) ■ = lim W W 

e^0+ ^^' 


^ c^(xg,x + e) - c'^(xg,x) af(x) 
|e| Cfe(x) ■ 


fc=l {x,e)(^El ' ' 

Since c®(xg, x + e) — c'^(xq, x) < |e|^^(x, e), we obtain 
[ Vxc(xo,-) -w^ < lim Y^ 

= / Otifk- 
Jn 
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Passing to the limit in 5 —>■ 0, we finally get 


/ Va;C(a;o,-) • < V] / afcCfe 

Jn Jn 

= / ^ ^Xki{x,Vkix)) 


dx. 


2.) First, using (14.91) with w = Ov for v G K'^) and an arbitrary scalar 

function 0 G C“(n, R), 0 > 0, we deduce from the homogeneity of z >->• z) 

that 

Va;c(xo,x) • ?;(x) < u(a::)), a.e. on fl. (4-10) 

Now let x be a Lebesgue point of both ^ and S/xc{xq,.), u G and take 
V G such that u = m in some neighbourhood of x. By integrating 

inequality (14.101) over Br{x), dividing by its measure and letting r —>■ O’*" we 
obtain 

Va;C(a;o,a:) • u < ^^{x,u), a.e. on fl. 

From iLemma 4.3l the desired result follows. □ 


4.2 The F-limsup inequality 

Given ^ G L^{0), we now prove the F-limsup inequality that is there exists a 
family G such that 

^ $ and limsup ) < J{0- (4.11) 

£—>0 + 

First show that for ^ continuous and then a density argument will allow us to 
treat the general case. 

Step 1 : The case where ^ is continuous. 

For every e > 0, {x, e) G , write 

We have 

—>■ IICIIp and /o(^®) -A loiO as e —>■ O'*". 

In particular, for e > 0 small enough, ||^®||e,p < 2||^||p and ^ in the weak 

sense of idefinition 3.11 We can proceed analogously to the construction (IT71) of 
c for the F-liminf. We define cf on the whole of fl x fl in a similar way and we 
also have the uniform convergence of c® to some c in (7(14 x fl) (passing up to 
a subsequence) and liminfg_,.o+ 4i(C'^) = J-Q^jicd'y so that to prove (j4.11|) it is 
sufficient to show that c > cj = cj. To justify this inequality it is enough to see 
that by construction for (x, y) G x one has 

c"{x,y)= inf f {a{t))dt > c^{x,y) 

o-eC|,, Jq 
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using the uniform convergence of (f to c we indeed obtain c > = c^. 

Step 2 : the general case where ^ is only L^{9). 

Let G C(0 X such that 

U - ?n|lp + I|C5„ - Cfiloo + \IoiCn) “ /o(^)| < “ 

and 

Un\\p<m\\p 

(existence is given bv ILemma 3.31) . For every n > 0,£ > 0 there exists G 
such that Cn Cn- Then there exists a nonincreasing sequence e„ > 0 
converging to 0 such that for every 0 < s < Sn we have 

\Ioien) - loiUl < mn) > hiU - - and < 2||e„||p. 

n n 

For e > 0, let rig = sup{n;e„ > s} and then we get —>■ C (IIC®l|e,p < 

2||^n||p < 4||^||p) as well as 

|/q(^'^) — Lo(C)I ^->■ 0 as e —^ O'*" 

ne 


and 

Itm > hiU) - - = f 

jQxn 

Since converges to c^, we then have 


1 

ne' 


liminf/i(^®) > /i(e) 


which completes the proof. 

5 Optimality conditions and continuous Wardrop 
equilibria 

Now we are interested in finding optimality conditions for the limit problem: 

inf J(^) := [ H{x,v,^{x,v)) d{dx,dv) — [ c^d'^, (5.1) 

46i+(S) JnxS'*-! JQxQ 

through some dual formulation that can be seen in terms of continuous Wardrop 
equilibria. More precisely, it is in some sense the continuous version of the 
discrete minimization problem subject to the mass conservation conditions (1^ - 
(IT^ . Write 

C = {{a,p) : a G Wi'°“([0, l],n),p G V, fl L°°([0,1])"^}, 


27 




where 


Per = < p : i e [0,1] H> p{t) e : ct(<) = ^ Vk{<T{t)) Pk{t) a.ei 
I fe=i 

We consider £ as a subset of (^([0,1],R'^) x £^([0,1])^ i.e. equipped with the 
product topology, that on (^([ 0 , 1 ],R'^) being the uniform topology and that on 
L^([0,1])-^ the weak topology. Slightly abusing notations, let us denote M\ {C) 
the set of Borel probability measures Q on (^([0,1],R'^) x £^([0,1])-^ such that 
Q{C) = 1. For cr G kF^’°°([0,1], n), let us denote by a the constant speed 
reparameterization of a belonging to kF^’°°([ 0 , 1 ], ff) i.e. for t G [ 0 , l],CT(t) = 
a{s~^{t)), where 

*(0 = TTT / l‘^(^)l with Z(ct) = / |(t(u)| du. 

Jo Jo 

Likewise for p G £°“([0,1])^, let p be the reparameterization of p i.e. 

Pkit) := I . Pfc(s~^(t)),Vt G [0,1], fc= l,...,fV. 

|cr(s 

We have p G Pa fl £°°([0,1])^ with ||p||li = HpHli. Define 

C := {{a, p) G C : \d\ is constant} = {{d,p), {a, p) G £}. 

Let Q G we define Q G A4^(£) as the push forward of Q through 

the map (cr, p) —>■ (cr,p). Then let us define the set of probability measures on 
generalized curves that are consistent with the transport plan 7 : 

2 ( 7 ) := {Q G M\(C) : (eo, ei)#Q = 7 }, (5.2) 

where cq and ei are evaluations at time 0 and 1 and (eo,e.i)^Q is the image 
measure of Q by (eo,ei). Thus Q G Q(j) means that 


(p(cr(0),cr(l)) dQ(a,p) := 


ip{x, y) d'y{x, y), \/(p G x 


Iflxfl 


This is the continuous analogue of the mass conservation condition (12.21) since 
Q plays the same role as the paths-flows in the discrete model. Let us now 
write the analogue of the arc flows induced by Q G 2 ( 7 ); for fc = 1,..., let 
us define the nonnegative measures on SI x by 


/_ ip(x,v) dm2(x,v) 
JcixS'^-^ 


p((7(t),Vk(cr{t)))pk {t)dt 


dQ{a, p), 


for every (p G (7(12 x Then the nonnegative measure on S2 x ^ 

t>e defined by 


/OxS"*- 


^drrfi = / L^(a, p) dQ{cr, p),y^ G C{n X 


d-l 


(5.3) 
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where for every (a, p) € £, 


N 

Ld<^,p) = 


^(cr(0, Vk{(T{t)))pk{t) dt= ^{cr{t)) ■ pit) dt, 


(5.4) 


with 

^i^it)) = • ■ •: d^id^VNicrit)))). 

Notice that L^ia,p) = Ldd,p) for every (cr, p) G £ and so m‘^ = rrfi 
for every Q G M.\iC). The p growth asumption p.7l) on H{x,v,-) can be 
reformulated by a g = p/{p — 1) growth on Gix,v,-). To be more precise, 
we will assume that gix,v, ) is continuous, positive and increasing in its last 
argument (so that G{x,v,-) is strictly convex) such that there exists a and b 
such that 0 < a < b and 


am^ ^ < g{x, v, m 

) < bimd ^ + 1) V (a:, u, to) G 14 X ^ x 1R+, 

(5.5) 

with q G (1, d/id — 1)). 

Then let us define 


2^7) : 

= {Q G Q(7) : rrfi G L«(f4 x S'^-\0)} 

(5.6) 

and assume 

2^(7) ^ 0- 

(5.7) 


This assumption is satisfied for instance when 7 is a discrete probability measure 
on n X n and q < d/(d — 1). Indeed, first for Q G AI)|_(TT^’°°([0,1], fl)), let us 
define ig G AI+(n) as follows 


f p diQ = f ( f ipia{t))\(7{t)\dt^ dQia) foi (fi G 

Jo. ./ivn°°([o,i],?2) \Jo J 


It follows from [4] that there exists Q G A4)|_(IT^’°“([0,1], fl)) such that (eo, ei)#Q 
7 and iq G L'^. For each curve a, let p'^ G Va such that Pkd) ^ C'|d(t)| (we 
have the existence thanks to [Assumption 4| ) . Then we write Q = (id, p')^Q. 

We obtain Q G so that we have proved the existence of such kind of 

measures. 

Let Q G Q'^i'j) and ^ and f be in (7(14 x S'^“^,R-|_), we have 


Lc(cr, p) - Lj(cr, p) dQ[a, p) = 


< 


iii^it)) -li^d))) ■ p{t) dt 


rrfi 9{dx, dv) 


dQicr, p) 


< lie - eilLp(e)ll"i'^l|L9(e)- 


So if ^ G L+(0) and dfin is a sequence in (7(14 x §‘^“^,]R+) that converges 
in LP{0) to ^ then is a Cauchy sequence in fi{£,Q) and its limit (that 
we continue to denote by L{) does not depend on the approximating sequence 
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(^n)ra- This suggests US to define in an L^{C,Q) sense for every ^ S T+(^) 
and Q S For every ^ S L^{&) and Q S Q'^(7), by proceeding as for 

Lemma 3.6 in [5], we have 

j ^ 9{dx,dv) = j L^{a, p) dQ{a, p), (5-8) 

dnxS'^-i Jc 

and 

c^(ct(0), tT(l)) < ^^(cr, p) for Q — a.e. {a,p)GC. (5.9) 

Hence using the fact that Q G Q‘^{'y) and (I5.8l) - (l5.9p . we obtain 

[_ _c^dj= [ c^{a{0),a{l)) dQ{a,p) < [ ^ ■ mP. (5.10) 

JnxQ. Jc dnxS‘i-1 

Let ^ G L^{0) and Q G it follows from Young’s inequality that 


/ H{x,v,^{x,v)) 9{dx,dv) 

JQxS'i-^ 

> I ^ ■ mP 9{dx,dv) — j G {x.v.mP {x,v)') 9{dx,dv) (5.11) 
so that we have 

inf J(^) > sup — / G [x,v,m^{x,v)) 9{dx,dv). (5.12) 

56i+(®) QeQ‘!(7) JnxS'i-^ 

The dual formulation of (EU then is 

sup — / G [x,v,mP {x,v)) 9{dx,dv). (5.13) 

Qea«(7) dnxS-^-i 

We can note the analogy between (j5.13ll and the discrete problem that consists in 
minimizing (12.41) subject to the mass conservation conditions (I2.2I) - (I2.3I) . Then 
we establish the following theorem, that specifies relations between (I5.13P and 
m and that gives the connection with Wardrop equilibria: 


Theorem 5.1. Under assumptions (15.51) and (EZl), we have: 

1. (15.131) admits solutions. 

2. Q G solves (15.131) if and only if 

p) dQ{a, P) = C4^(o-(0), o-(l)) dQ{a, p) (5.14) 

where fqix, v) = g v, m^{x, v)^ . 

3. Equality holds : inf dSTj = sup (15.131) . Moreover if Q solves (15.131) then 
f-Q solves (j5.ip . 
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It is the main result of this section. To prove it, we need some lemmas. 
First, let us start with a preliminary lemma on the Vk's that is a consequence 
of [Assumption 4| 

Lemma 5.1. For all subset I C {1,..., N}, it is all or nothing, that is, we are 
in one of the two following cases : 

1. 0 € Conv{{vi{x)}i^i) for every x € fl, 

2. 0 ^ Conv{{vi{x)}i^j) for every x € fl. 

Moreover, there exists a constant 0 < ^ < 1 such that for all subset I C 
{!,...,A^} that is in the second case, there exists Ux £ Conv{{vi{x)}i^i for 
all X € such that 

II 

vAx) ■ > 5 for all i G I. 

\Ux\ 

Proof. We will use the fact that SI is connected. The first property is obviously 
closed since the Vk’s are continuous. Let us now show that the second one is 
closed. Let I C {l,...,Af}, assume by contradiction that there exists a sequence 

{a;r!,}rt>o £ ^ converging to a; £ SI such that 0 ^ Cn = Conv({ui(a;n)}iG/) for 
every n > 0 and 0 £ C = Conv({tii(a;)}ig/). So there exists {Ai}ig/ such that 
= 0, Ai > 0 and ~ Without loss of generality, we can 

assume that the A^’s are positive. Then we have that Kvi{xn) ^ 0 

and converges to 0 as n —>■ +oo. Let fdn > 0 such that \/3nv'^\ = 1 then /3„ 
converges to +oo. Thanks to [Assumption 4[ with ^ = 0 if 

i £ /, 1 otherwise, there exists such that |zf| < C and 

E^e7 i(xn) = PnVi{xn) for all 7 £ J and n > 0. Then we obtain that 
— zf)vi{xn) = 0. But for n large enough, we have that /3„Ai — z" > 0 
for every i G I, which is a contradiction. 

For a subset / C {1,... ,iV} that is in the second case, Conv({ui(a:)}ig/) is 
contained in a salient (pointed) cone for all cc £ fl. For all x £ SI, we can think 
of Ux as being in the medial axis of this cone. Since the are continuous and 
SI is compact, we have the desired result. □ 

Now let us define the following sets 

f ^ 

= I {<^,P) & ■ '^Pkit) < C\a{t)\ a.e. t G [0,1] 

I k=l 

for some constants C > 0. Now let us notice that we can simplify the problem 
([5. 131) with the following lemma: 

Lemma 5.2. For a well-chosen constant C > 1, one has 

inf f G {x,v,nfi{x,v)) 9{dx,dv) 

QeQiG) InxS‘i-1 

= inf If G (x,v,ni'^ (x,v)) 6{dx,dv) : Q (c^ ^ = 1 

QeQ9(7) ^ 
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Proof. We set C' = 1/6 where 5 is given bv ILemma 5.11 Let (cr, p) S C. We 
will prove that there exists p £ Va such that for all t G [0,1], Pfc(t) < Pk(i) 
for all fc = and — C"|(T(t)|. Let t G [0,1] such that 

J2k=iPk{t) > C"|tT(t)|. Let us denote I the subset of {l,...,iV} such that 
for every k G I,pk{t) > 0. First, if 0 G Cone({rifc(cr(t))}feg/), there exists a 
conical combination of 0 

^ = 0 
k^I 


with the Afe’s > 0. Then we write Pk(t) = Pk{t) — XXk (we take Xk 
where 


A = min 

kGl:Xk^0 



0 for fc ^ I) 


We set I the subset of / such that for every k G I,Pf^{t) > 0. We restart with p 
and we continue until 0 ^ Conv({z;fc(CT(t))}j,gj). Let u be as in ILemma 5.11 for 
I = I and X = cr(t). Then we have 


N 


N 


\&{t)\ > a{t) • w = ^ Pk{t)vk{a{t)) pf,{t) 


k^l 


k^l 


so that Pki^) ^ C"|(T(f)|. For Q G we denote by Q G 

the push forward of Q through the map (cr, p) i—^ (c, p). Then we have 771*5 < . 

Since G{m,v,-) is nondecreasing, we have : 


/ ^ 



0{dx,dv) < 




G (a:, u, 777‘5 (x, w)) 9{dx,dv). 


□ 

To prove that the problem (15.131) has solutions, a natural idea would be to 
take a maximizing sequence {Qn}n>o for (15.131) and to show that it converges 
to Q G that solves (15.1311 . For this, we would like to use Prokhorov’s the¬ 

orem which would allow us to obtain the tightness of {Qn} and *-weak conver¬ 
gence of {Qn} to a measure in Unfortunately, the space G([0,1], x 

L^([0,1])^ is not a Polish space for the considered topology (because of the weak 
topology of L^{[0, 1])). So we will work with Young’s measures in order to apply 
Prokhorov’s theorem. Let us define the set C = G([0, IjjR^^) x tPi(R'^ x [0, Ij) 
where for a Polish space {E,d), we set 

^i{E) = G M.\{E) : J d{x,x') dp{x) < -boo for some x' G E 

We equip C with the product topology, that on G([0,1],]R‘^) being the uniform 
topology and fPi(]R‘^ x [0,1]) being endowed with the 1-Wasserstein distance 

TUi(p^,p^) := min|y d(a;i,a;2) dp{xi,X2) ■ p G n(p\p^) 
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where E = x [0,1], d is the usual distance on E, G ^i{E)'^ and 

, fi'^) is the set of transport plans between and /r^, that is, the set of 
probability measures ^ on £ 1 ^, having and y? as marginals: 


f ExE 


(p{x)dii{x,y) = / and / ip{y)dfi{x,y) = / ip{x)dy 

Je Jexe Je 


’ix), 


(5.15) 

for every (p G C{E,'B.). The set C is a Polish space (see [T]). Let us denote by 
A the Lebesgue measure on [0,1] and let us consider the subset 5 of C : 

S = {{a, Et ^ X) : a G 1T^’°°([0,1], 11), r-t (g) A G ^i{E), vt G 911^ a.e. t} , 

where for t G [0,1] and cr G 1T1’“([0, 1],11), 


mi = 


Vt G 


M\{ 


N 

) : supp ut (Z \ \ R+Vk{cr{t)) and d{t) = v dvt{v) 

k=i "" 


The Young measures vt®X are the analogue of the decompositions p GVa- 
For the general theory of the Young measures, see for instance m 

Let us define the set of probability measures on curves (u, vt ® A) that are 
consistent with the transport plan 7 : 

Y( 7 ) := {X G M^(5) : (eo, ei)#X = 7 }. (5.16) 

This is the analogue of (15.21) . Let us now write the analogue of (given by 
dssj) as follows: 


/ ^di^ = [ L^{a,K)dX{a,K)y^GC{nxS‘^-\R) (5.17) 

where for every (cr, k = X) G L, 

L^{a,K)=J (^J^^^(^a{t),-E^\v\dEt{v)^ dt. (5.18) 

Then let us define 

X«( 7 ) := {X G Y( 7 ) : G L«(fl x (5.19) 

Let X G By the same reasoning as for Q G Q'^( 7 ), if ^ G L+(0), we 

denote by the limit of the Cauchy sequence in L^{S, X) for any sequence 
(^ra)ra converging in LP{9) to We may write the analogue of the problem (I5.13P 

sup — f G (^x,v,i^{x,v)) 9{dx,dv). (5.20) 

Lemma 5.3. One has sup (|5.13l) = sup (|5.20ll . 
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Proof. Let Q G Q'^(j) and a G VL^’°“([0,1], ri). For p = (pi,..., p^) G Va, we 
define the Young’s measure t'f 0 A as follows : 




Pkit) ^ 


where |p(t)|i = Pk{t) for every t G [0,1]. We consider the measure G 

^9(7) defined by 

I p dX^ = I p{a,v^ ® \) dQ{a, p), for all € (7(5, R). 

Js Jc 

Since we have rrfi = we immediately get sup (15.1311 < sup (15.201) . 

For the converse inequality, let X G we build G Q‘^(j). Let 

(a, 0 A) G S, recalling that one has supp t't C U^i ^+'‘’k((^(i)) for t G [0,1], 

we define p" G as follows 


Pfe(0 = [ 1^1 di^t{v), for all fc = 1,..., 


and p'' = {pi,... ,p((f) if the Ufc(cr(t))’s are pairwise distinct. Otherwise, let us 
decompose {1,..., TV} = Uj=i dj where the IkS are pairwise disjoint and such 
that for all j = 1,. .., s and k G Ij,Vk(<j(t)) = Vj where the Uj’s are pairwise 
distinct. Then for all j = 1,..., s and k G Ij, we set 


pm 



|u| dvt{v). 


The element p'^ is in Va- Similarly, we set 


J pdQ^— J p{a, p'^) dQ{a,h't ® ioT all (/? G (7(£, R). 

From the fact that it follows that sup (15.131) > sup (15.201) . □ 


Let us notice that with the previous proof for (cr, 0 A) G 5, we may build 
7 as a sum of Dirac measures : 




N 


pm s: 


where p^ is given in the previous proof. Therefore it follows from the same 
reasoning as in the proof of ILemma 5.21 that we may take p'^ G Va such that 
for t G [0,1], ^ C"|d-(t)|. Moreover, we can choose (cr, p) only in £ 

with \a\ constant. Then the new measure \p''{t)\i ^{\pit)\i‘Vki<^{t))} that we 

continue to denote by 7* by abuse of notations is in SOl^. Let us define 


= {(cr, r-t 0 A) G 5 : supp Vt n R+Ufe(cr(t)) = {p{t)vk{cr{t))} 

with p(t) < C'\&{t)\ for fc = 1,...,TV and t G [0,1]} 
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and 


5 = {(cr, ( 8 ) A) G 5 : \&\ is constant}. 

For X G Ad}_(5), we denote by X G the push forward of X through 

the map (cr, vt ® X) ^ {a, vt ® A). Then we have . Since G{m, v, •) is 

nondecreasing, we may consider only the measures X G 0 5) for the 

problem (15.201) . 

We now adapt the proof in [2] . In particular we have to generalize Lemmas 
2.7 and 2.8 in [ 8 ], this becomes 

Lemma 5.4. For every ip G C(n x L^p is l.s.c. on S for the topology 

defined above. 

Proof. Let (cr, z/j 0 A) G 5 and (cr", z/" 0 A) be a sequence converging to (cr, z/* ® 
A) G S. Then by definition, we have 

(g) A) = ^ ip |^cr"(t), Ini di^fiv)^ dt. 

^ a in (7([0,1]) so that </j(cr"(-), |^) converges strongly in L°°. Since vf 
narrowly converges to vt and the function (t,n) i—>■ p{a(t), ■i}(|-)|n| is the upper 
limit of (t,n) i—>■ (/ 5 (cr(t),-pji-) min(|n|,n), that is continuous and bounded, as 
n —>■ + 00 , we obtain the desired result. □ 

Lemma 5.5. Let (X„)„ G Ad}_(C')^ he such that Xn{S'^ ) = 1 for every n and 
there exists a constant M > 0 such that 

sup / l{(j) dXn{cr, zcj (g) A) < M. 

n JS 

Then the sequence {Xn)n is tight and admits a subsequence that weakly-k con¬ 
verges to a probability measure X such that X(S) = 1. 

Proof. For every K >0, let us define the following subset of 5^ 

Sk = |(cr, zcj (g) A) G 5 : |cr| < if and supp vt < Bc'if | 

where C is the constant given bv ILemma 5.21 Let us show that Sk is relatively 
compact in S. First, the set {cr G IF^’°“([0,1], fl) : cr if-Lipschitz continuous} 
is compact in C([0, l],n) thanks to Ascoli’s theorem. The set of probability 
measures with support in Bc'k is compact due to the Banach-Alaoglu-Bourbaki 
theorem. Let a sequence (cr”, z/” (g) A) G (5ic)^ converging to (cr, vt ® X) G 5, 
prove that (cr, zc* (g) A) G Sk- 

1.) supp Vt C \Jk=i^-i-'Okicr(t)). 
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First let us notice that the function (p : {x,v) !->■ dist(u, M+Ufe(x)) is 
continuous on x We then have : 

[ ([ ‘fiv,a{t)) diytiv)] dt = f (f p{v,a{t)) diytiv)] dt 
Jo \JB<i / Jo xJBaiK J 

p{v,cr^{t)) dt 


^C'K 

r-1 


= lim 

n—>-+oo 
= 0 . 


0 \J Br 


So p{x,v) = 0 dvt ® dt-a.e. and the support of t't is in IJ^j^ K.+Ufc(CT(t)) for 

tG [0,1], 

2-) d{t) = ^^^vdvt{v). 

By definition, for n > 0 and (s, t) G [0,1]^, we have : 


a"(t)-a"(s) = 


V dv^{v) ® \ = (wlsc 


(g) A). 


Obviously, the sequence {CT”(t) — cr”(s)}„>o converges to a{t) — cr(s) (since cr” 
uniformly converges to a). For the term in the right-hand side, it is sufficient to 
take a sequence {$e}£>o in C'6(]R'^ x [0,1])'*^ converging to (v, t) !->■ (g)l[s_t] 

in X [0,1]) as £ —>■ O’*’. 

Now let us justify the tightness of (X„)„: 





< ^|((T, t-t (g) A) G 5 n 5*^ 

-I- Xn r-t (g) A) G 5 n 5*^ : supp(r't) ^ Sc'_ftr|) 

< 2X„ ({(a,:.* 0 A) G 505^' : |d| > x}) 


< 2X„ ({(ff, (g) A) G 5 : 1(a) > K}) 

< ^ y l{<x) dXn{a, Vt <g> A) 

^ 2-— —y 0 as K —y -t-oo. 


Due to Prokhorov’s theorem we can then assume that passing up to a sub¬ 
sequence, (Xn)n weakly-* converges to Af G M\{C). It remains to show that 
Af(5) = 1. For iG > 0, let us define the closed set 


Sk = {(a, r'i (g) A) G 5 : l{a) < K and supp vt C Bck} ■ 

It follows from the previous computation, the fact that the measures Af„ are 
concentrated on S and Portmanteau’s theorem that 

I = limsupAf„(5) < limsup Af„(5/f) -I- limsup Af„(5\5K) 

n n n 

<X{Sk) + ^. 
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Letting K tend to oo, we then obtain X[S) = sup^X(5/f) = 1. □ 

Lemma 5.6. Let (Xn)n be a sequence in that weakly star converges to 

some X S If there exists i € + x sueh that i^”- weakly-* 

converges to i in Al + (f2 x then we have i^ < i. 

The proof is similar to that of Lemma 2.9 in [5]. 

Proof. (of ITheorem 5.11) 

Let us prove the existence of solutions for the problem (15.131) . Thanks to 
ILemma 01 we consider the problem (15.201) . Due to (15.51) the value of problem 
(I5.20|) is finite. Let (X„)„ be a maximizing sequence of (|5.20|) . Since i^ > i^, 
we can assume Xn = Xn for all n. Still from (15.51) it follows that is 

bounded in L'^{9). So, passing up to a subsequence, we can assume that {i^'')n 
weakly converges in L'^{6) to some i. Moreover, since is bounded in L‘^{9) 

so in L^{9), we have 


sup / l{a) dXn{a,i^t <Si X) < sup { { Ivldi'tiv)] dt\ dXn{(J,vt <Xi X) 

n Js n Js \Jo \Jr'^ / / 


= sup 


di^" < + 00 . 


n JQxS'^-^ 


Since Xn = Xn, we can deduce from ILemma 5.51 that, up to a subsequence, 
[Xn)n weakly-* converges to some X G Af?|_(5). Using the fact that X{'-/) is 
weakly closed, we see that X G and ILemma 5. 61 then imply that < z so 

that X G Since G(x,v,-) is convex and nondecreasing, we then have 


/OxS'^- 


G{x,v,i^(x,v) 9{dx,dv) < / G{x,v,i{x,v)) 9{dx,dv) 

1 ./nxS‘i-1 

G{x, V, z^" {x, v)) 9{dx, dv), 


< liminf / 


which proves that X solves (|5.20l) . Thus as mentioned in the proof of ILemma 5.31 
there exists Q G such that nfi = i^ and so Q is a solution of (I5.13|) . 

The reasoning for the last two statements is similar to that of Theorem 5.1 

in [2]. □ 


A natural question is to investigate the discrete problems corresponding to 

dH i.e. 

e m^(x,e)' 


mf lel'G X,-, 




(5.21) 


subject to the mass conservation conditions (I2.2p - (j2.3l) and convergence of prob¬ 
lems (I5.2ip in some sense to the continuous problem 


inf f G{x,v,nrfi{x,v)) 9{dx,dv). (5.22) 

QeQ(7) dnxS'^-i 
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Let m® = {in^{x,e))(x,e)eE'^ and w® = {w^{u))a^c^ solve the discrete prob¬ 
lem (15.211) . Let (j = {xq, ... ,Xi(o.)) G C® (identified with the piecewise affine 
curve defined on [0,L(cr)]). For every fc = 0,... ,L{a) — 1, let us denote by ik 
the integer such that {xk,Xk+i — Xk) G . Then let us define p’^ G L°°([0,1])^ 
where for all t G [fc, fc -|- 1[, 


Pi (t) = 


\a{k -b 1) - cr(fc)| if i = ik, 
0 otherwise. 


We will define a discrete measure over where 


= {{a,pn--aGC^}. 


Write as follows 

QE ._^d/2-l ^ W%a)da0p,^ 

o-eCs 

as well as 

Q" ^ w'^{a)S^^pa 

o-ec= 

where tr G 1], fl) is the constant speed reparameterization of the path 

a. Notice that for every ^ G C{fl x we have L^{a,p'^) = L^{a,p’^) 

so that . Let us also observe that the measure contains all the 

information on (m^,w®). 

Especially for the following theorem, we make a stronger assumption. 

Assumption 10. There exists a funetion C : M+ i—?> R.^ sueh that C{e) —>■ 1 as 
e —>■ 0'*' and for every e > 0, {x, e) G E^,C(s)s < \e\ < e. 

In particular, this hypothesis is satisfied in our three classical examples since 
arc length is constant for £ > 0 fixed. 

Theorem 5.2^ Under the previous assumptions, defining Q® as above, up to a 
subsequenee, {Q^)s > 0 weakly eonverges to some solution Q G Q'^(j) of (I5.22|) 
in the sense that 


/ ^i<^,p)dQ‘^{(T,p) ^ / ^{a,p)dQ{a,p), 

/C([0.1],R‘i)xLi([0,l])" dc([0.1],R‘i)xLi([0,l])" 


as £ 


0+ for every $ G ^(^([0,1], x Li([0,1])^, R). 

Proof. By duality, from ITheorem 5. H and Corollary 3.1[ it follows that the value 
of (15.2111 converges to that of (15.2211 and in particular, due to the q growth 
condition (j5.5|l on G{x,v,-), m® is bounded for the discrete norm. In the 
same manner that in the proof of [Corollary 3.1| and Section 4.1 we can see that 
there is some m G L+ such that (a:, e) —» ^ weakly converges to m in L'^ 


nuuu uiiau —r |g|d/ 

in the sense of idefinition 3.1 1 (up to replacing p by q) and 


/nxS'i- 


G{x,v,m{x,v)) 9{dx,dv) <\iiRini |e|‘^G|a;, ■ 


(x,e)£E^ 


'{x,e) 


j|d/2 


(5.23) 
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Let ^ S C(f2 X S'^ ^,R+), recalling (I2.3L (13.41) and (15.3p . rearranging terms, we 
have 


IQxS‘‘~ 


^(x,v)drrfi {x,v) = / L^{a, p) {a, p) 


= e-^/ 2 -i 




= e-^/ 2 -i 


^fc+i 

II IZ 

cr^C^ k—0 

L{a)-1 

H w"(a) ^ 

<7^C^ k—0 

L(a )-1 

(?(o'(^) 


+ 1) - o-(fc)|dt 




a{k + 1 ) — (T(fc) 


o-eC« 


fe =0 


[<T(fe).<T(fc+i)] V’k(^ + i)-c^(^)l 

a{k + 1 ) — cr(fc) 


+ C>(wj(£:)) 


a{k + 1 ) — o-(fc)| 

^£,d/2-i ^ ^2; ®'j I ^ |e|w^(cr) 

(a;,e)G£^^ ^ ' \ /rcr^e-F-r- /t 


+ 0(u’5(f)) |(t(/ 5 + 1) - a(A:) 




k cr€C^:[x,tc+e]Ccr 




e \ m^(x, e) , ,, 

’^+ 0 (zc 5 (e)) 


;|!i /2 


where is a modulus of continuity of From [Assumption 10| and the fact 
that {x, e) —>■ weakly converges in L'^ to m in the sense of idefinition 3.11 

it follows that weakly star converges to m. Arguing as previously, we find 
Q S Mi{C) such that, up to a subsequence, {QOe weakly converges to Q and 
mP < m. We easily have Q G ■ indeed, for every p G C'(K'^ x K^,K), we 

have 


/ dQ{a,p) = lim / p{a{0),a{l)) dQ^a, p) 

Ic Jc 


= ^hm ^ u;-(a)^M0),a(l)) 


o-gC' 


(x.y)GAr '2 

= lim V p{x,y)-fPx,y) 

E-> 0 + 

{x,y)£N^^ 


p dj. 


/QxQ 
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Using (15.231) and the fact that G{x,v,-) is nondecreasing, we get 




G{x, V, nrfi{x, v)) 9{dx, dv) < 


/ G{x, V, m{x, v)) 6{dx, dv) 


< liminf |e|‘^G 

e->0+ , '^7^ 

{x,e)GE^ 


f e m^{x,e)\ 


Since the right-hand side is the value of the infimum in (15.221) . we obtain the 
desired result. □ 


6 The long-term variant 


Instead of taking the transport plan 7 *^ as given in the discrete problem, we now 
consider the case where only its marginals are fixed. More precisely, there is 
a distribution of sources ft = sinks /+ = X^xgW' 

which are discrete measures with same total mass on the set of nodes (that 
we can assume to be 1 as a normalization) 

E E /+(?/) = 1- 

xGN^ yGN^ 

The numbers fl{x) and /+(a;) are nonnegative for every x € N^. 

With the same notations as in the short-term problem, we have almost the 
same definition of an equilibrium as in [definition 2 .ll we must change the mass 
conservation condition ( 12 . 21 ) as follows 

f-ix) := E /+(y) •= E 

for every {x,y) G x N^, where G^ . (respectively G^y) is the set of loop-free 
paths starting at the origin x (respectively stopping at the terminal point y). 
Moreover, the transport plan now is an unknown. Similar arguments apply to 
this case, the equilibrium is a minimizer of the functional defined by (12.41) but 
now subject to (ED) and (IQ) . We shall then state the analogue of the dual 
formulation (12.51) 


inf 


E H‘^{x,e,t'^{x,e)) 

{x,e)GE^ 


inf 

7'en(/i,/') 


E l’'{^^y)T!e{x,y) 


{x,y)^N^'^ 


( 6 . 2 ) 

where n(/i, /^) is the set of discrete transport plans between ft and /^, that 
is, the set of nonnegative numbers ( 7 ^( 2 ;, 2/))(a; such that 


E 7''(2:,y) = E 7 ^( 2 :, y) = /+(y), v(x,y) x Nt 


We assume that the hypotheses made in Subsection 3.1 are still satisfied, except 
that we replace [Assumption 7|by 


40 











Assumption 11. /I and weakly star converge to some probability measures 
f- and /+ on SI; 


lim ^ {(p{x)fl{x)+tp{x)fl{x)) = f >pdf_+[ tjjdf+y{(p,tp) £ C{n)'^. 

e-o+ Jn Jn 

Writing as in (13.81) . we can now reformulate (16.21) 


inf (6.3) 

eeK#® 


where /o(C^) defined by (I3.10|) and 


) := inf 

7=en(/i,/=; 


E 

{x,y)GN^'^ 




|d/2 




(6.4) 


' ( 2 ,e)Co’ 


It is an optimal transport problem. The limit functional then reads as the 
following variant of (13.161) 

.^^(6 := do(6 - where Fi(^) := inf / G (6.5) 

7en(/_./+)7nxn 

As previously, /q is defined by (13.121) and by (13.151) . !!(/_,/+) is the set of 
transport plans between /_ and /+ (see (15.151) 1. We then have the following 
F-convergence result : 

Theorem 6.1. Under the same assumptions except \As sumption ''/\ replaced by 
\Assumption 11\ the family of functionals F® defined by (16.31) T-converges (for 
the weak -topology) to the functional F defined by (16.51) . 

Same arguments as for Theorem 6.1 in [5] apply here. 

In the same manner as in ISection 5[ we can see that the problem (16.51) has 
a dual formulation that is 


sup 

Qea<iU-J+) 


IcixS"^- 


G{x, V, m'^{x, v))0{dx, dv), 


( 6 . 6 ) 


where 


Q^(/_, /+) := {Q G Mt(£) : eo#Q = f-,ei^Q = /+, G T«(0)} 

= U 

7en(/_,/+) 

If we assume that /+) ^ 0 and that (15.51) is still true, one can reformulate 

[Theorem 5.11 for the long-term models as follows : 

Theorem 6.2. We have : 
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1. Problem (EID admits solutions, 

2. Q G S'*(/_,/+) solves (16.61) if and only if 

dQ{a,p) = J^c^^{a{0),a{l)) dQ{a,p) 

where £;q{x,v) = g{x,v,nrfi{x,v)) and moreover, 7 := (eo,ei)^Q is a 
solution of the optimal transport problem: 

, /_ _CiQ{x,y)d-/{x,y). 

7Gn(/_,/+) jQxn 

3. There is no duality gap : the infimum of dSlD equals the supremum of 
(16.6p and moreover, if Q solves (16.6|) then solves (16.5|) . 

Problem (16.61) is studied in [S]. It is showed that problem (16.61) is equiva¬ 
lent to another problem that is the variational formulation of an anisotropic, 
degenerate and elliptic PDE : 

J — div (f\/Q*{x,'Vu{x))) = f in n, 

VG*{x,Vu{x)) ■ vq = 0 on dGl, 

with Q* being a function. In particular, if the function g in (13.61) is of the 
form g{x, Vk(x), m) = ak(x)m'^~^ + Sk for every x G Gl,k = 1,..., N and m > 0 
where the constants Sk are positive and the weights Ok are regular and positive, 
then we have 


G*ix, z) = _ SkCk{x)Y, for every a; € fl, z € 

fe=i ^ 

where bk = {akCk)~'^■ This case is interesting since numerical simulations can 
be performed as shown in [9]. 
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